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Introduction.—Tidal forces, and the rotations, cause the fluid components 
of a binary star to deviate from a spherical shape. Over the surfaces of 
such distorted stars, gravity varies. The shape and surface gravity of each 
component depend upon its size, the rate of rotation of the matter within it 
and its distribution of density, as well as upon the separation and masses. 
Chandrasekhar * * 4 has given formulae by which both the shape and 
surface gravity can be computed when the distribution of density is poly- 
tropic. Elsewhere,’ the present author has given formulae by which the 
shape may be computed corresponding to any distribution of density, and 
these formulae agree exactly with Chandrasekhar’s, as they should, when 
the distribution is polytropic. In the present note, formulae are given for 
the surface gravity corresponding to any distribution of density. The sur- 
face gravity is of interest because the bolometric surface brightness is be- 
lieved to be proportional to it.® 

We use the notation and procedure of reference 5, in which formulae are 
given (on page 460) involving the second, third and fourth order harmonic 
distortions, for the shape of a rotationally and tidally distorted fluid star 
having any distribution of density. The method there followed employs a 
type of analysis first developed by Laplace. Denote the masses of the 
stars by m, and m2; their centers of mass by Q; and Q2; the distance QQ, 
by R (which, in an elliptical orbit, is a function of the time). The total 
potential at a point A near star 1 is 


V=U+); (1) 


where U is the gravitational potential of the distorted star 1, and where V, 
is the disturbing potential: 
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Va=V,+ V,. (2) 


Here 
mG — ( r\" 
y= pears 
, R »» ( i) P,,(cos 6) (3) 
is the tidal potential due to star 2, and V, is a rotational term 
Y= - 3 on'r*Px(cos 6’) (4) 


arising from the rotation of the matter in star 1 with a uniform angular 
velocity w, about an axis normal to the plane of the orbit. In these ex- 
pressions 7 is the distance QA, G is Newton’s constant of gravitation, P, 
is the Legendre coefficient of order n, 0 is the angle AQ,Q2, and 6’ is the co- 
latitude. It is shown in reference 5 that if one ignores small quantities of 
the order of the squares and cross products of the Y’s, later to be defined, 
then corresponding to each term of the form 


rr (B) 
in'the disturbing potential, where n > 2, there is a contribution to the super- 
ficial distortion given by 


Qn+1 a*t!} 


n + 1n(a) mG ”" (6) 





Y, =c 


and a contribution to the external gravitational potential of star 1 given by 


i 2 1 — n,(a) 
n ttl n + n,(a) n° 





(6) 


Here a is the mean radius of star 1, and the external shape of star 1 is 
given by 
r=a(1+ DY,). (7) 


Different P,,’s need not be referred to the same axes; and the Y, and U, 
always have the character of the surface harmonic, P,, present in the dis- 
turbing term (B) that causes them. 

In these formulae 7,,(a) is a number that depends upon m and upon the 
distribution of density in star 1; that is to say, upon the stellar model. It 
may be computed, for any particular stellar model, by solving the first order 
differential equation 

dn, 


6 
roe t mt — 14 — mn + 1) + “lm + 1) = 0 (8) 





P= 





VoL. 27, 1941 ASTRONOMY: T. E. STERNE 95 


subject to the boundary condition that at r = 0, », = m — 2; and then 
n,(@) is the value of n, at r = a. To define p and p,, precisely it is really 
necessary to consider the surfaces of constant density inside the distorted 
star. A surface of any given density, p, is characterized by a parameter 7, 
its mean radius. In the differential equation p is the density at the distance 
r from the center, and p,, is the mean density interior to 7, in a spherical 
distribution of matter, not necessarily in equilibrium, in which the density 
is the same function of 7 that p, in the actual distorted star, is of 7. The 
specification of the stellar model may be regarded in practice as defining the 
p in equation (8), and thus p,,, as a function of r. For any particular 
model, equation (8) can always be solved easily by numerical integration so 
as to yield the number »,,(a); for certain models, analytical solutions can be 
found. 
It follows that the rotational distortion of star 1 is given by 











a ) am (9) 
a 
and the tidal distortion by 
pret n+ % (10) 
where 
273 
Yior = — Ae = Py(cos 0’); Y, = — A,s"t4P.(cos 6); (11) 
3m,G my, 
and where 
2n +1 1 
= R); A, = ——; PF: = —(3x? — 1); 
£= (0/R); A, = 5; Pals) = 532" — 1) 


Pi) a =(5x3 < S: Payee =(35e4 — 30x + 3). (12) 


Values of As, As and A, are given by Chandrasekhar® ‘ for the polytropes; 
he obtained the values not by solving the present equation (8) but by 
special analytical procedures. It has been mentioned that the squares and 
cross products of the Y’s have been neglected in the present treatment. 
The smallest square or cross-product term to be neglected is the square of a 
Y;, term, and this square is of the same order of magnitude as Y;. For 
most purposes the present degree of accuracy, which is the same as Chand- 
rasekhar’s for polytropes, is sufficient. If a still higher degree of accuracy 
were desired, it would be necessary to take into account not only the 
squares and cross products, and harmonic distortions with orders higher 
than the fourth, but also corrections to V, arising from the distortion of 
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star 2 produced by its own rotation and the tidal action of star 1. In 
particular, there is little to be gained by considering the square of a Y2 term 
without considering Y;, and little to be gained by considering Y; without 
taking into account the squares of the Y2 terms. To the order of accuracy 
of the present treatment, the various Y’s and U’s are additive. 

Consider, now, a term in the disturbing potential 


cr" P,,(cos 6) (B) 


where 6 is measured from any appropriate axis—in the case of a tidal term 
it must be measured from the line of centers, and in the case of a rotational 
term from the axis of rotation. Let us consider the effect of this term 
upon the surface gravity. Considering only this single disturbing term 
we shall have for the total potential in the external vicinity of star 1 
v= ood + U,, + cr"P,,(cos @) 

in which the first term on the right is the potential due to the undisturbed 
star 1. The corresponding value of apparent gravity, g, at the surface of 
star 1 is of course 


dv 


dn 


evaluated at the surface, where dz is measured in the direction of the out- 
ward normal. The value of gravity being continuous at the surface, we do 
not need the value of W at interior points; and g is the resultant of the dis- 
turbing force and the attraction of the distorted star 1. The outward nor- 
mal makes a small angle with the radius vector r. It is easy to show that 
this angle is of the order of magnitude of Y2 or smaller, so that its cosine can 
differ at most from unity by a quantity of the order of magnitude of the 


square of Y2. To the present accuracy, therefore, 
_dv 
dr 


evaluated at r = a(1+ Y,,). Making use of the expressions already given 
for U, and Y,, one finds after some reduction that 


— = —[1 + »,(@)] Y, (13) 
0 
where 
5 
a= (14) 
a 


and is a constant. 





= He 


= - = 
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More generally, if a number of terms of the form (B) in the disturbing 
function are to be considered, then 


cB = - Dll + nC) ¥,. (15) 

This equation has a simple meaning: for each separate harmonic distor- 
tion there is a corresponding variation of surface gravity, such that the 
fractional variation of surface gravity is proportional over the star’s surface 
to the fractional variation in the radius vector. The factor of proportion- 
ality is 1 + 7,(a) and depends on the stellar model and on the order of the 
distortion. 

A Theorem.—From equations (5) and (13), and from the form of the 
term (B), it follows that the apparent surface gravity, corresponding to 
any harmonic distortion of order , obeys the law 





aa -(2 + Ve + (n — 1)¥y, (15’) 
So n 

where ¢, is the ratio of the outward normal component of the disturbing 
force at the surface to gp. Equation (15’), not used here, enables one to 
compute the value of apparent gravity over the surface of any tidally and 
rotationally distorted fluid mass from a knowledge of its shape and of the 
values of the disturbing forces at the surface, even when the internal con- 
stitution is unknown. If the actual surface distortion is analyzed into 
zonal (or tesseral) harmonics of different orders, Y,, then the total (g — 
£0) /go is merely the sum of the right-hand member of (15’) over all the zonal 
(or tesseral) harmonics that are present. When 1 is set equal to two, 
equation (15’) becomes Clairaut’s well-known relation that is used in find- 
ing the earth’s ellipticity from measures of gravity. I am indebted to Dr. 
Russell for having pointed out to me that some general relation—such as 
the exceedingly general equation (15’)—was necessarily contained in the 
analysis of this Note and in that of reference 5. 

Surface Gravity on a Component of a Binary Star.—In the practical appli- 
cation of equation (15), there are two Y’s of the second order, namely Y,,, 
and the tidal Y2. There are also the tidal Y3 and Y, of the third and fourth 
orders, respectively. Thus, in a binary star, the surface gravity g of com- 
ponent 1 is given by 


aids 4 
as = —[1 + m(a)] Yee — 2{1 + 24(0)] ¥, (16) 
0 2 
where the Y’s are given by equation (11). 
The 7,,(a)’s are related to the k,,’s, of importance in the theory of apsidal 
motion, in the manner indicated in reference 5. For all models, 
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e+ 1 — in (2) 
eae 
The “‘K”’ of Russell’s paper* is the present ke; Russell does not consider the 
third and fourth order harmonic distortions. The quantity 


4, [1 + 1,(@)] = 2 + 2 — 2k,(n — 1) 


gives the dependence of the range of variation of surface gravity, over the 
surface of a star, upon the distribution of density and upon the order of the 
harmonic distortion. The 7,(a)’s can be found for any density distribution 
by solving equation (8), as has been mentioned; and they can be found 
without further work for polytropes by using Chandrasekhar’s A’s and the 
relations (12). For the reader’s convenience we list in the accompanying 
table the values of k, A, 1 + n(a) and A[l1 + n(@)] for polytropes. The 
polytrope of index zero is a homogeneous model; that of index five is a 
completely concentrated, Roche, model. The values for any model, 
polytropic or not, should lie between the values corresponding to the poly- 
tropic indices zero and five, since these are limiting distributions of density. 
Observational evidence, obtained principally by Cowling,’ Sterne® * and 





POLYTROPIC INDEX 0 1 1.5 2 3 4 5 
ke 3/, 0.2599 0.1446 0.0741 0.0144 0.00134 0 
As 5/, 1.5198 1.2892 1.1482 1.0289 1.00267 1 
1 + m(a) 1 2.2898 2.8784 3.3546 3.8596 3.9867 4 
As[1 + m2(a) ] 5/, 3.4802 3.7108 3.8518 3.9711 3.9973 4 
ks 3/g 0.1064 0.0540 0.0244 0.00368 0.00024 0 
As 7/7, 1.2129 1.1079 1.0488 1.00736 1.00047 1 
1 + ns(a) 2 3.7713 4.3183 4.6743 4.9489 4.9967 5 
As[1+7:(a) | ™/, 4.5742 4.7842 4.9024 4.9853 4.9991 5 
ke 1/, 0.0602 0.0281 0.0116 0.00140 0.00007 0O 
A 3/, 1.1205 1.0562 1.0231 1.00281 1.00014 1 
1 +m(a) 3 5.0821 5.5211 5.7968 5.9748 5.9987 6 
A.{1 + m(a)] 9/, 5.6385 5.8314 5.9307 5.9916 5.9996 6 


Russell® from the apsidal motions, suggests that actual stars are consider- 
ably concentrated. They appear to be sufficiently concentrated for values 
of A, 1 + »(a) and A[1 + n(a)] to be taken from the last, Roche, column of 
the table with sufficient accuracy for most practical purposes. In many 
practical applications, moreover, it will be legitimate, in the absence of 
conflicting observational evidence, to set the constant angular velocities w, 
and w: equal to the mean orbital motion 27/P where P is the period. 
Summary.—in Monthly Notices, 99, 460 (1939), the author gave formu- 
lae for predicting the shape of a rotationally and tidally distorted fluid 
component of a binary star. The component could have any distribution 
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of density, and harmonic distortions through the fourth order were taken 
into account. In the present note formulae, valid for any distribution of 
density, are given for the surface gravity. A tabular compilation is made 
of the numerical coefficients that are needed to compute the apsidal mo- 
tions, distortions and surface gravities of some polytropic models. 


1 Chandrasekhar, M. N., 93, 390 (1933). 

2 Chandrasekhar, Jbid., 93, 449 (1933). 

3 Chandrasekhar, Jbid., 93, 462 (1933). 

4 Chandrasekhar, Ibid., 93, 539 (1933). 

5 Sterne, Ibid., 99, 451 (1939). 

6 For references see Russell, Ap. J., 90, 641 (1939). 
7 Cowling, M. N., 98, 734 (1938). 

8 Sterne, Ibid., 99, 662 (1939). 


NOTES ON BINARY STARS. II. ON THE LIGHT OF ECLIPSING 
BINARIES OUT OF ECLIPSE 


By THEODORE EUGENE STERNE 


Harvarp COLLEGE OBSERVATORY 


Communicated January 8, 1941 


Introduction and the Second Order Harmonic.—In Note I in these Pro- 
CEEDINGS,! formulae were given for finding the surface gravity of a rota- 
tionally and tidally distorted component, of a binary star, having any 
distribution of density. Here we give a brief derivation of the effects of 
the distortion, and of the variation over the surface of the surface bright- 
ness, upon the observable light of the components outside of eclipse. The 
method used is simple and general, and will have applications to other phe- 
nomena involving the light of a distorted star having a non-uniform surface 
brightness. 

A theorem was proved by von Zeipel* and Chandrasekhar* that over a 
tidally and rotationally distorted star, the bolometric surface brightness 
must be proportional to gravity if the star is in hydrostatic equilibrium. 
The argument can be put in a simpler and more physical form than the 
proofs given by the two previously mentioned authors. The equation of 
hydrostatic equilibrium is dP = pd¥, where P is the pressure, p is the 
density and W is the total potential asin Note I. It follows that P and p 
must be functions of WV alone in a star that is in hydrostatic equilibrium 
under the influence of its own gravity and of any disturbing potential. 
Hence the temperature, 7, must also be a function of Y alone and so must 
be the opacity since that is a function of T and p. The outward flow of 
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heat is given by the product of the temperature gradient and a function of 
the temperature, density and opacity. The function is a constant over 
any equipotential surface in the star, and thus the outward flow of heat 
across such a surface is proportional, over the surface, to the temperature 
gradient. Since the equipotential surfaces are isothermal surfaces, the 
temperature gradient is proportional over such a surface to the gradient of 
W and the latter gradient is merely the value of apparent gravity, g. Thus 
the outward flow of heat is proportional, over any equipotential surface 
within the star, to gravity; and if one takes for the selected equipotential 
surface the outer surface of the star, or a surface close within it, it follows 
that the bolometric surface brightness is proportional to apparent surface 
gravity. 

The preceding physical demonstration reveals the limitations of the 
theorem. The theorem need not hold if there are convection currents near 
the outer surface, for then the equation of hydrostatic equilibrium is in- 
applicable. The equation will be in error by amounts depending upon the 
acceleration of the stellar matter near the selected equipotential surface. 
It does not seem likely that this source of error can be very serious; a more 
serious one is that an appreciable transport of heat by material currents 
may render the flow of heat no longer proportional, over the selected equi- 
potential surface, to the temperature gradient. It is difficult to predict 
theoretically, for actual stars, how much the bolometric surface brightness 
deviates from the strict proportionality, to surface gravity, that would hold 
under conditions of perfect hydrostatic equilibrium. Moreover, even if 
exact proportionality holds for the bolometric surface brightness, it can 
hardly hold for the surface brightness in a single wave-length, for the 
quality of the radiation must vary over the distorted surface, with the sur- 
face brightness. 

Let us take for the variation of surface brightness, as it is influenced by 
the presence of a harmonic distortion of order n, the expression 


J= J, — £¥a)- 


Here J is the surface brightness as influenced by this distortion, J, is the 
average surface brightness and Y,, (defined in Note I) is such that the radius 
r is given by 

r=1,(1 + Y,) (17) 


where 7, is the average radius. If the surface brightness is bolometric, and 
if the theorem is applicable, then the coefficient 8, is exactly equal to the 
quantity 1 + »,(a) whose value las been found, as a function of the stellar 
model and of the order m of the harmonic distortion, in Note I. To take 
account of possible deviations from the theorem, the safest course appears 
to be to set 
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B, =y[1+n,(@)] (18) 


where y is the coefficient by which the full variation, predicted by the 
theorem, must be multiplied to obtain the actual variation of surface 
brightness. If the theorem is applicable, y = 1; if the surface brightness is 
independent of surface gravity, y = 0. The coefficient y is the same as 
Russell’s,4 and some observational evidence that the bolometric surface 
brightness follows the law y = 1 was presented by Zdenék Kopal to the 
New York Academy of Sciences in a lecture given on October 20, 1939, and 
is contained in a paper, now in press, that he has been kind enough to 
show me. 

We also take into account the darkening of an apparent stellar disc 
towards its limb by using the usual law of limb-darkening that renders the 
apparent surface brightness proportional to 1 — x + x cos y, where y is 
the angle between the outward normal and the line of sight. Thus 


J = Jl — 8,¥,)(l — & + % cos 4), (19) 


where J, is the average normal surface brightness, for each harmonic com- 
ponent of the distortion separately. 

As in Note I, we shall ignore the squares and cross products of the Y’s, 
and to this order of accuracy the fractional effects of the different distor- 
tions are additive. As pointed out previously,’ it is useless to consider 
small quantities of the order of the square of the oblateness or higher, un- 
less harmonic distortions of orders higher than four are considered; here 
we shall consider harmonic distortions of orders up to and including the 
fourth. To the present accuracy, no squares or cross products of the Y’s 
can appear in any results. The rotational distortion Y,,, hardly con- 
cerns us, for it can introduce no observable effects in the light out of eclipse. 
Out of eclipse, the aspect of the rotational ellipsoid does not change since 
the direction of the axis of rotation does not change. Thus the observed 
variation in light does not involve any effect proportional purely to Yyo. 
Moreover, Y,,.¢ can combine with another distortion Y,, only in the form of 
a product Y,., Y,, and we have agreed to ignore such small cross-product 
terms. Since we are not concerned with the light during eclipse, we can 
thus legitimately ignore the rotational distortion so far as changes of light 
are concerned. By similar reasoning we see that we can consider the effect 
of each tidal distortion Y,, separately, as though it alone were present, and 
we proceed first to consider Y2, the most important tidal distortion. 

The tidal distortion Y2 causes the surface of star 1 to become approxi- 
mately a prolate ellipsoid, with its long axis along the line of centers.* 
The radius to any point on the surface is given by 
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Y2 = (2/3)eP2(cos 5’) 
so 


r = 1 + (2/3)eP2(cos 5’) (20) 


correctly to the first order in the tidal oblateness, ¢, if the mean radius is 
taken as unity. Here 6’ and y are ‘‘astrocentric’”’ spherical polar coérdi- 
nates in which 6’ is the angle between r and the line of centers. The polar 
coérdinates of the outward normal to the surface at the given point may be 
denoted by 6 and y; and then in terms of the “‘astrographic”’ codrdinates 6 
and y it is easy to find expressions for the radius vector, 7, and the element 
of surface, dS. Since 


£. = —2e sin 6’ cos 6’ (21) 
rdé 
one has to the first order in e 
6’ = 6 — 2esin 6 cos 6 (22) 
and so r = 1 + (2/3)eP2(cos 5). Moreover, since the outward normal 


makes an angle to the radius vector that is proportional to ¢, the cosine of 
the angle is unity to the order of «. Thus 


dS = r* sin 6'di’dy 
= [1 — (8/3)eP2(cos 5)] sin ddédy 
or 
dS = [1 — (8/3)eP2(cos 8)|dQ (23) 


correctly through terms of order e. Here dQ is an element of solid angle in 
the normal system 4, y. 

Take a right-handed system of axes X, Y, Z such that the Z-axis is the 
axis of the orbit and the Y-axis is in the direction of the ascending node of 
the relative orbit of star 1. Let the inclination be denoted as usual by 12, 
and let 6 be the true orbit longitude of star 1 measured from superior con- 
junction. When the orbit is circular, this @ is the usual 6 of eclipsing star 
theory. Take also a set of rotating axes X’, Y’, Z’ with Z’ towards the 
earth, and X’ in the plane of the sky and towards the apparent position of 
star 2. The angle a between the Z’-axis and the line of centers is given by 


cos a = sin 7 cos @. (24) 


We employ a set of spherical polar angles y, ¢ in the primed system, such 
that y is the polar angle measured from the Z’-axis. In the primed sys- 
tem, the direction cosines of the direction y, ¢ are sin y cos ¢, sin y sin ¢, 
cos y; those of the line of centers are sin a, 0, cos a. Hence 


cos 6 = sin a sin y cos g + cos a cos y. (25) 
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The angle y is just the angle involved in the limb-darkening; and dQ = 
sin ydydg. Accordingly, the light received at the earth from the element 
of surface corresponding to dy and d¢ on star 1 is proportional to 


dl = J,{1 — x + xcos y] [1 — (2/3)e(4 + Be)Pe(cos 8)] 
‘sin y cos y dydg_ (26) 


and the total light is proportional to this quantity integrated over all 
values of y from 0 to 7/2 and of g from 0 to 2x. The integration is easily 
effected since cos 6 is given by (25); the result is that the light is propor- 
tional to 





, 9 Tee 
Jf 1 — 5 — Zeb + Pu(cos «(+7 * + Ax) 


so that one can write for the light of star 1 


eee Ld 1 _ (Ba + 4)(15 + 2) 
30(3 — x) 





P,(cos 2) | (27) 


correctly to terms of order e where L, is a constant, and if we are concerned 
only with the changes in light caused by the changing geometry} during an 
orbital period we can write this as 


(Be + 4)(15 + x) 
20(3 — x) 





L= Ly 1-. sin? 7 cos? | (27’) 
where L,’ is a new constant. Equation (27’) is equivalent to Russell’s 
equation (14) of reference 4, since his term that involves the rotational 
ellipticity is a mere constant whose effect can be absorbed, to his accuracy 
and ours, in the constant L,’. 

If ‘‘gravity darkening”’ is ignored and there is no darkening towards the 
limb, then $2 and x are both zero and (27’) agrees to the order of € with 
Russell’s equation (41) of reference 9 in which his ‘‘e’”’ is an eccentricity and 
his ‘‘e?”’ is our 2e. If gravity darkening is ignored and there is complete 
darkening towards the limb then f; is zero and x is unity and equation (27’) 
agrees as it should with Russell and Shapley’s equation in the text on page 
400 of reference 8. In accordance with the definition of 8. and the results 
in the previous Note, the full gravity effect (y = 1) corresponds to a value 
of 62 equal to four for the Roche model and to a value (which in any par- 
ticular case can be exactly computed as in Note I) of nearly four for any 
strongly concentrated model. As noticed by Russell,‘ therefore, a full 
gravity effect must in practice nearly double the variation in light, observed 
outside of eclipse, that arises from the principal tidal distortion Y2. 

An equation similar to (27) holds, of course, for star 2; and the observed 
light outside of eclipse is the combined light. 
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The Third and Fourth Order Harmonics.—For the third order harmonic, 
one has 
r = ¢3P;(cos 6’) (28) 


where ¢; is a numerical coefficient that may be found as described in Note I. 
By the same procedure as before, one finds that 


dS = [1 — 10e;P3(cos 5)|dQ (29) 
and thus the light of star 1 outside of eclipse, as influenced by the third 
order tidal deformation, is given by 

2x 2/2 

SS (l — x +x cos y| [1 — (10 + 83)Ps(cos 4)| sin y cos ydyd¢ 

0 0 
and thus the light of component 1 varies in accordance with the equation 
(10 + Bs)x 
4(3 — x) 








L= Lt — 6 P;(sin 7 cos 0) | (30) 


Here 63 = y[(1 + m:(a)]. 
For the fourth order harmonic tidal deformation, 


= 1 + «&P,(cos 6’) (31) 
where «, is a numerical coefficient given in Note I, and one readily finds that 
dS = [1 — 18e,Ps(cos 6)]dQ. (32) 


Thus the light of star 1 outside of eclipse, as influenced by the fourth order 
tidal distortion, is proportional to 
2x 4/2 


SS (Wl — x +x cos yJ| [1 — «(18 + Bs)Pa(cos 6)] sin y cos ydyde 
0 0 
so that the light varies in accordance with the equation 


L=L,{i + os Ba 2. Bi i i cos 6). (33) 
8(3 apie? x) 





We notice that no variation of light arises from the third harmonic unless 
there is some limb-darkening, nor any from the fourth unless there is in- 
complete limb-darkening. These particular properties of the odd and even 
harmonics higher than the second were noticed by Russell in a paper” on 
the light-variation of asteroids. 

If there is a complete gravity effect, then from Note I we see that 63 and 
8, are nearly five and six, respectively, for strongly concentrated stars. 
The exact values can be computed for any particular stellar model. The 
full effect of the variation of gravity is thus in general to multiply the varia- 
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tions of light outside of eclipse, arising from the third and fourth order dis- 
tortions, by factors of 3/2 and 4/3, respectively. 

If the right-hand members of equations (27), (30) and (33) are written as 
L,(1 + Az), Lo(1 + As) and L,(1 + d4), then the combined variations give 
for the light of star 1 the quantity Z,(1 + A» + As + Ay). A similar result 
holds for star 2. 

One must remember, when combining the lights of stars 1 and 2, that the 
“@” for star 2 differs by 180° from the @ for star 1. Thus the second and 
fourth order variations of star 2 are given by equations (27) and (33) with 
appropriate coefficients, but in the third order variation of star 2, given by 
equation (30) with appropriate coefficients, P; must be replaced by —P3 if 
for 0 one employs the @ of this Note, which is the 6 of star 1. If the com- 
ponents are identical, one sees that no odd-order harmonics can appear in 
the combined light whatever may be the law of darkening; for the effects 
of the odd harmonics of star | are then precisely neutralized in the combined 
light by those of star 2. On the other hand, the even-order harmonics of 
stars 1 and 2 must always reinforce each other. 

The effects of reflection must of course be added to the effects considered 
here, in order to obtain a complete and accurate representation of the light 
outside of eclipse. 

Systematic Effects —The quantity L, is the light of star 1 free from the 
effects of the tidal distortions. Since the integral of a surface harmonic 
over a sphere is zero if the order differs from zero, L, is thus dependent 
merely upon the luminosity of star 1 and its distance from the earth—ex- 
cept for reflection and the small effect considered shortly. The quantity L,’ 
of equation (27’) contains, however, a systematic tidal effect arising from 
the constant term in P2(cos a); and it should not be used without correction 
to obtain a luminosity free from tidal effects. It can easily by corrected 
by the equation 


(34) 





a Ly 1 _ (Be + 4)(15 + |. 
60(3 — x) 

The rotational distortion does not influence the variation in light outside 
of eclipse to the present order of accuracy, but it tends to introduce a sys- 
tematic error in the luminosity. This error can be allowed for by the rela- 
tion (obtained by the same general procedure as equation (27)) 


(B2 + 4)(15 + x) 
30(3 — x) 





L, = [1 + €rot P(cos | (35) 


where L* is free from the effects of the unequal distribution of brightness 
over the surface of the rotational ellipsoid. Here ¢,,, is the rotational 








106 ASTRONOMY: T. E. STERNE Proc. N. A. S. 


oblateness, equal to (—3/2) times the coefficient of P2 in Y,., given in 
equation (11) of Note I. 

Summary.—Equations are obtained for the light, as influenced by rota- 
tional and tidal distortions, of a component of an eclipsing star outside of 
eclipse. The results are valid for any distribution of density, and all har- 
monic distortions are taken into account with orders up to and including 
the fourth. The conditions of applicability of von Zeipel’s theorem about 
surface brightness are examined, and some systematic corrections are con- 
sidered to the observed absolute magnitudes of rotating stars. 


* In the limit where the orbital period of a binary star is long compared with the free 
periods of spherical harmonic oscillation of the components, the tidal deformations are 
given by the equilibrium theory of the tides. See Cowling, reference 5, and Sterne 
reference 6. For the most complete discussion of the ratio of periods see Sterne, refer- 
ence 7, pages 664 and 665. In the present Notes, the equilibrium results are used in 
which the line of centers is a common axis of all the tidal harmonics. 

+ Changes in light caused by changes in e during the motion will be studied in a future 
Note. 

1 Sterne, Proc. Nat. Acad. Sct., 27, 93 (1941). 

2 yon Zeipel, M. N., 84, 702 (1924). 

3’ Chandrasekhar, Jbid., 93, 589 (1933). 

* Russell, Ap. Jour., 90, 641 (1939). 

5 Cowling, M. N., 98, 734 (1938). 

6 Sterne, [bid., 99, 451 (1939). 

7 Sterne, Jbid., 99, 662 (1939). 

8 Russell and Shapley, Ap. Jour., 36, 404 (1912). 

9 Russell, Jbid., 36, 54 (1912). 

10 Russell, Jbid., 24, 1 (1906). 


NOTES ON BINARY STARS. III. THE ELLIPTICITY EFFECT 
IN ECCENTRIC ECLIPSING BINARIES 
By THEODORE EUGENE STERNE 
HARVARD COLLEGE OBSERVATORY 
Communicated January 11, 1941 


In Note II we derived? equations for the variation in the light of a binary 
star outside of eclipse, as influenced by the distortion and by the non- 
uniform distribution of surface brightness. The tidal distortion and the 
distribution of surface brightness depend in part, as has been shown in 
Notes I' and II, upon the separation. When we computed the light in 
Note II, we considered the distortion and surface brightness to stay con- 
stant during the motion, as they would in a circular orbit. If the orbit is 
eccentric, however, the separation changes and causes the distortion and 
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distribution of surface brightness to change, in the period of the orbit. 
Here we give the simple equations for the effect of these changes upon the 
light outside of eclipse.* 

The light of star 1, as influenced by the harmonic tidal distortion Y2 of 
order two, is given outside of eclipse by equation (27) of Note II. In that 
equation, «¢ is the tidal oblateness of the second order; it is three-halves of 
the coefficient of P: in the expression (11) for Y2in Note I. Let v denote 
the true anomaly, othe longitude of periastron, and e the eccentricity of the 
relative orbit, of star 1. Then if the mean distance between centers of 
mass is A, (A/R) = (1 + ecosv)/(1 — e*). But by equation (11) of Note 
I, « varies as (1/R)*; thus « = e’(A/R)* where ¢’ is the value of ¢ at the in- 
stants of mean separation and is given by 








e’ = (3/2)(me2/m) Ae(a/A)?*. (36) 
Thus by equation (27) of Note II one has for the light of star 1 
(Bet 4(15+x), . 2. (1 + ecosv)* 5} 
L= Ig E —~¢'- 60(3 — x) (3 sin? 7 cos? @ — WT a (37) 
where 


6=o+7-— 90° 


and where Ly is the light of star 1 free from the effects of the second order 
tidal distortion. 

In equation (37), e’ is a constant during the motion; and so is the whole 
coefficient e’(82 + 4)(15 + x)/60(3 — x), whose value may be denoted by 
S.,. The combined lights of stars 1 and 2 outside of eclipse, as influenced 
by the changing tidal distortion of the second order, will thus obey the 
relation 


(1 + ecosv)? 


38 
er, (38) 


Leombinea = 1 — S2(3 sin? 7 cos? @ — 1) 
where S2 is a mean of S,, and S.» (for star 2), weighted with the Ly’s of the 
two stars, and should be observationally determinable. In (38) we have 
taken the combined Ly’s as unity. 

In the same way, the third and fourth order harmonic tidal distortions 
produce light variations for star 1 of 


, 4 
L=I, E — S;,(5 sin*z cos’ 6 — 3 sini cos 6) dof eceea| (39) 
—e 


and 


: pee pare 5 
i= I E + S,,(35 sin*7 cost @ — 30 sin? 7 cos? 6 + 3) eee 
—e 


(40) 
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respectively, outside of eclipse. Here the coefficients are given by the 
equations, obtained from Notes I and II, 


= (10+ Be _ ms ) (18 + 44)(1 — 2) 
ele a (5)! "8(3 — x)’ ™ a (4 64(3 — x) 


(41) 





and remain constant during the motion. Equations (39) and (40) hold for 
the combined lights of the two stars, if the Zy’s are set equal to unity and the 
S's are replaced by S; and Sy, the weighted means, respectively, of S;, and 
— S32, and of S,; and S,». 

In this Note, we content ourselves with giving the theory of the “‘distor- 
tion’’ terms only, in the light. To the changes in the combined light of the 
components due to the distortions, there must, of course, be added the 
changes in light due to the reflection effect, if one wishes to obtain an exact 
and complete representation of the observed light outside of eclipse; and 
in the reflection effect, the changing separation must be taken into account. 

We notice that equations (38) and (40) for the combined lights both 
predict small differences (in the same sense) between the combined lights 
at the ascending and descending nodes (9 = 270° and 6 = 90°), unless@ 
is either 90° or 270°. The reflection effect, being greater near periastron 
than near apastron, should always reinforce this difference, which the third 
tidal harmonic does not influence. 

The author has been informed that Dr. Kopal has a paper in press on the 
subject of this Note—the effect of changes in figure upon the light outside 
of eclipse. Dr. Kopal’s investigations, and the author’s, have been inde- 
pendent. The present Note contains, substantially, the contents of a 
paper that the author read on December 30, 1940, before a meeting of the 
American Astronomical Society in Philadelphia. 

Summary.—The effects are predicted of changes in figure, occurring 
during a revolution in an eccentric orbit, upon the light outside of eclipse. 
Ordinary and ‘gravitational’ darkening are taken into account, and the 
results are applicable to any distribution of density. 


* Cowling® has considered the effect of the changes in the second order tidal distortion, 
and the author‘ has considered the effect of the changes in the second, third and fourth 
order tidal distortions, upon the longitude of periastron. The theory of the effects of the 
changing distortions upon the light is, needless to say, considerably simpler than the 
theory of their effects upon the orbit. 

1 Sterne, Proc. Nat. Acad. Sci., 27, 93-99 (1941). 

2 Sterne, Ibid., 27, 99-106 (1941). 

3 Cowling, M. N., 98, 734 (1938). 

‘Sterne, M. N., 99, 451 (1939). 
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BIVALENT STRUCTURE IN THE FLY, MELOPHAGUS OVINUS 
L. (PUPIPARA, HIPPOBOSCIDAE) 


By KENNETH W. COOPER 


DEPARTMENT OF BIOLOGY, PRINCETON UNIVERSITY 


Communicated December 20, 1940 


The customary absence of genetic crossing-over in male Drosophila made 
difficult the acceptance of any hypothesis of meiosis in which the chiasma is 
the condition sime qua non of association of homolegous chromosomes in 
bivalents at the first meiotic metaphase. Darlington (1931) recognized 
the difficulty inherent in Drosophila males and suggested that meiosis of the 
male brachycerous Diptera is exceptional not in the absence of chiasmata, 
and hence of crossing-over, but in that the chiasmata present are invariably 
double, reciprocal,! and located in genetically inert regions of the chromo- 
somes. However, this view was later retained only for the sex chromo- 
somes, for his own work (1934) as well as that of Keuneke (1924), Metz 
(1926) and Bauer (1931) seemed to show that the association of homologous 
autosomes in higher Diptera is derived from an exaggerated expression of 
the forces of somatic pairing. Since physical connections have not been 
generally observed between the four chromatids in autosomal bivalents, 
Darlington holds the forces of repulsion between homologous chromo- 
somes to be in exact equilibrium with forces of attraction at meiotic meta- 
phase. The absence of genetic crossing over in the male brachycerous fly, 
therefore, is currently believed by Darlington to be due to an anomalous 
meiosis in which the autosomes are held together as bivalents only by an 
abnormally intense expression of the attractive forces of homologous genes, 
and the sex chromosomes by the occurrence of nearly invariably reciprocal 
chiasmata in the genetically inert but homologous regions. 

Of the suborder Brachycera of the Diptera, especial interest lies in the 
section Pupipara.? Not only have the Pupipara remained unexamined 
cytologically, but they represent perhaps the biologically most specialized 
derivatives of the Schizophora—the very series to which Drosophila be- 
longs. The present paper gives a preliminary account of certain critical 
features in the meiosis of the male of Melophagus ovinus L.—the sheep 
tick, or ked—with a brief discussion of their significance for the general 
theory of the meiosis of male Brachycera. It will be shown that not only 
does Melophagus possess a unique type of autosomal bivalent for flies, 
but that the structure of these autosomal bivalents sheds new light upon 
the interpretation of the structure of sex chromosome bivalents in Dro- 
sophila. 

The gonial metaphase plate of the male Melophagus ovinus (Fig. 2) shows 
18 chromosomes.* The diploid set may be resolved into a group of eight 
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large autosomes, eight small autosomes and an unequal XY pair. The 
large autosomes are all of similar size and each bears a trabant-like arm 
proximal to its kinetochore. The small autosomes likewise are all of similar 
proportions, and each bears a trabant-like region distally and, perhaps, 
another proximal to its kinetochore. Whereas the autosomes all show a 
more or less paired association on the mitotic metaphase plate, the sex 
chromosomes evince no such tendency (Fig. 2). 

The most significant peculiarity of the somatic mitoses (including gonial 
and embryonic cells, and follicle cells of the ovary) is the derivation of the 
paired condition of mefaphase from a proximal association or pairing in the 
large autosomes at prophase (Figs. 1, 4). The paired regions are terminal,‘ 
proximal to the kinetochore and may show faintly staining wisps connect- 
ing the homologues.’ Other than by proximal terminal pairing, no ten- 
dency toward somatic pairing is shown by these prophase chromosomes. 
The expressions of affinity are localized in the proximal trabant-like re- 
gions, and the remainders of homologous chromosomes appear mutually 
indifferent. 

In the small autosomes the condition is not so clear. Frequently both 
ends of homologous small chromosomes appear paired, while the mid- 
regions bow away from one another as though in a state of active repulsion. 
Not uncommonly small autosomes are paired at one end only; in such 
cases the paired ends are clearly those distal to the kinetochore. At no 
stage do the sex chromosomes evince any tendency to undergo somatic 
pairing in the male (Figs. 1, 2, 4). 

In the prophase of the first spermatocyte there has been found no inter- 
pretable evidence of the early synaptic stages of meiosis. In this respect 
Melophagus resembles other cyclorrhaphous Brachycera.’ At diakinesis, 
however, the situation is strikingly different from that described for other 
flies such as Drosophila which have exceedingly small chromosomes. Both 
classes of autosomes clearly show chromatic (feulgen-staining) connec- 
tions as the mode of bivalent association. Although the large autosomes 
condense in a common karyosomic cluster and only infrequently show their 
structure clearly, the small autosome bivalents are quite evenly distributed 
throughout the nucleus and easily studied. At Jate diakinesis the small bi- 
valents are held together by distal connections, whereas in the large auto- 
somes there are similar proximal unions (Fig. 5). As the first metaphase 
plate forms, the repulsions of the kinetochores become especially evident 
and clearly locate these points of association (Fig. 6). Before the plate is 
fully formed the associative strands of the bivalents become attenuated 
(Fig. 6) and may finally appear to have disjoined. The result is a meta- 
phase plate in which the bivalents appear to be held together only by Dar- 
lington’s forces. 

The sex chromosomes during all stages of meiosis remain entirely apart 
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from one another and show no evidences of association (Figs. 3, 5, 6). 
Their behavior is unique among Brachycera, but resembles that of the 
supernumerary chromosomes of 7ipula (Bauer, 1931) and the sex chromo- 
somes of Fungivora, Mycetophila and Exechia (Frolova, 1929), for they re- 
main in the equator during early anaphase (Fig. 7). Only after the auto 
somes have attained the poles do the sex chromosomes reductionally sepa- 
rate. 


5 N " 6 "se 


SPERMATOGENESIS OF Melophagus ovinus L.; FiGuRES SEMI- 
DIAGRAMMATIC 


Figure 1, gonial prophase. Figure 2, gonial metaphase. Figure 
3, polar view of first meiotic metaphase plate. Figure 4, gonial 
prophase. Figure 5, diakinesis; one small autosome bivalent 
missing. Figure 6, prometaphase of first meiotic division; chro- 
mosomes have been spread across equator for purposes of clarity. 
Figure 7, anaphase of first meiotic division. It must be noted 
that figures 1 and 4 are approximately one-third again the magni- 
fication of the remaining figures. 


The second spermatocyte division is interesting in that the small auto- 
some univalents separate first, then the sex chromosomes divide equation- 
ally and lastly the large autosomes undergo anaphase. Although there is 
thus a definite order in which each chromosome class separates, the chromo- 
somes all reach the poles simultaneously. 

The chief importance of the striking chromosome behavior of male 
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Melophagus lies in four points. First, forces of intimate somatic pairing 
appear rigorously localized in their expression to relatively minute termi- 
nal regions of the autosomes in mitosis (Figs. 1, 4). Second, these special 
regions provide the mechanism of diakinetic association of homologous 
autosomes (Fig. 5). Third, the diakinetic connections between autosomal 
homologues, like chiasmata, give origin to the metaphase associations of 
the first spermatocyte (Fig. 6). Fourth, the sex chromosomes have evolved 
a mechanism of segregation employing neither chiasmata nor somatic 
pairing (Figs. 3, 5, 6, 7). In this respect the sex chromosomes behave 
similarly to those of certain Nematocera. 

Clearly, as in Darlington’s view for Drosophila and other Brachycera, 
the pairing of the autosomes at meiosis in Melophagus must probably be 
considered to be due to exaggerated somatic pairing. It differs from the 
reported association of Drosophila in that in Melophagus the pairing is 
strictly localized terminally, and results in a pronounced, continuous bridge 
between homologous autosomes of a bivalent. Contrary to the case of 
Drosophila, repulsions of the kinetochores at meiosis in Melophagus take 
place much prior to metaphase and the paired connections between homolo- 
gous autosomes become greatly drawn out. As the metaphase association 
of Melophagus follows one in which these obvious physical connections exist 
between the homologous autosomes, it is accordingly not unreasonable 
that, despite an optical void between the members of a bivalent, the asso- 
ciation is still maintained by their sharing a portion of common pellicle or 
by connection of homologous halves to the same spindle elements. The 
anomaly of repulsion of homologues in exact equilibrium with persistent 
specific attraction among all four chromatids of a bivalent need not be in- 
voked to account for the configuration of these bivalents of Melophagus 
at metaphase.’ 

The most significant features in the meiosis of Melophagus with regard 
to the general theory of meiosis in flies, however, lies in the following three 
points. (1) The homologous autosomes are only in contact over a short 
portion of their length, and this is assuredly not the whole of their homolo- 
gous regions (Figs. 5, 6). (2) Parts of homologues not paired—in spite 
of their homology—iie apart at diakinesis (Fig. 5). (3) At prometaphase 
(Fig. 6) and metaphase there appears to be a state of tension between the 
spindle attachments and the points of association.* Briefly, the autosomal 
bivalents of Melophagus give evidence of the very characteristics of struc- 
ture and behavior that Darlington holds to be indicative of chiasmata in 
the sex chromosomes of Drosophila. Nevertheless, in Melophagus they are 
probably all simulacra of chiasmata. As the relatively minute pairing 
regions of Melophagus seem to be merely chromosomal foci to which so- 
matic pairing forces are restricted, and which by exaggeration of these forces 
hold bivalents together at meiosis, we are not as yet forced to assume with 
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Darlington a cytological need to suppose reciprocal chiasmata in the inert 
regions of the sex chromosomes of Drosophila.’ Interstitial localization of 
somatic pairing forces at meiosis would give the figures of association that 
Darlington, Dobzhansky and others have published for the sex chromo- 
somes of Drosophila, without involving chiasmata. 


I wish to express my gratitude to Prof. Franz Schrader for suggesting 
this problem to me, to Prof. Schrader and Prof. Th. Dobzhansky for ex- 
tended critical discussion of the results and to Mr. J. Monahan, foreman of 
the New York Stock Yards, for making it possible for me to easily obtain 
large numbers of Melophagus. 


1 Darlington (1931, see dia. III) first supposed the kinetochore to be included between 
the reciprocal chiasmata. This view can no longer be maintained, and it must now be 
supposed that the reciprocal chiasmata of the sex chromosomes are both on one side 
of the kinetochore. 

2 Darlington (1934) followed Imms, who relegated the Pupipara to a separate series of 
Brachycera. 

3 This is contrary to Lassman’s report (1936) that the maturating egg and somatic 
cells show the chromosome number of Melophagus ovinus to be 16 (n = 8). Gonial 
plates of both male and female Melophagus show the diploid number to be 18. Ap- 
parently Melophagus shares with Anthrax sinuosa Wied. (Metz, 1916) the largest diploid 
number yet reported in the Diptera. 

4 In figures 1 and 4 the chromosome pairs at 12 o’clock and 11 o’clock, respectively, 
appear to be in lengthwise association. Actually the members of each pair are diver- 
gent. 

5 These connecting threads strongly recall those figured by Bauer (1931) for mitotic 
(gonial) prophases of Tipula. 

6 It must be asserted with Metz (1926), Huettner (1930) and Dobzhansky (1934) that 
it is not at all certain that the customary synaptic stages have been conclusively shown 
to be absent in Brachycera. Indeed, if Darlington is correct in his identification of 
chiasmata in the sex chromosomes of Drosophila, these chromosomes at least should 
reveal the customary sequence of stages. 

7 Presumably Darlington (1934, especially Fig. 36 of the side view of a diakinetic 
autosomal bivalent; 1937, especially pages 372, 498) has observed terminal connections 
or contacts between autosomal homologues at one extremity of the diakinetic bivalent 
of Drosophila. Whether these close associations, believed by Darlington to be due to 
pairing of the kinetochores, play a réle similar to the autosomal pairing segments in 
Melophagus can only be conjectured, but it should be pointed out that Stevens (1908, 
Figs. 67, 69, 70, 71) and Dobzhansky (1934, Figs. 12, 18, 23, 35) have also published 
figures which call to mind the phenomena here described for Melophagus. It is interest- 
ing to note that Darlington (1931) used Stevens’ figures in support of his earlier con- 
tention that fly autosomes form bivalents by reciprocal chiasmata. Each of these 
types of connection between autosomal bivalents deserves reinvestigation. 

’ 8 The wording is nearly identical with that of Darlington (1934) describing the sex 
chromosome bivalents of Drosophila. 

® It must be pointed out that Prof. Th. Dobzhansky believes that new evidence he has 
discovered in Drosophila duncani Sturt. supports Darlington’s reciprocal chiasma hy- 
pothesis for the sex chromosomes of Drosophila. 
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SOME EFFECTS OF THE GENE W’ (DOMINANT SPOTTING) IN 
MICE 


By ELIzABETH FEKETE, C. C. LITTLE AND A. M. CLOUDMAN 
Roscog B. JACKSON MEMORIAL LABORATORY 
Communicated January 3, 1941 


The gene W (dominant spotting) was identified by Little (1915) as hav- 
ing lethal effect when present in a homozygous condition. It was reported 
by de Aberle (1927) that among the offspring of matings in which both 
parents were heterozygous for the W gene, an inherited anemia occurred in 
certain animals. This appeared about the 16th day of intra-uterine life, 
became progressively more severe and resulted in the death of the animals 
within 10 days after birth. Little and Cloudman (1937) described a new 
dominant spotting mutation (W’) which probably represents a new allele 
of W. They stated that the homozygous offspring of this mutation may 
have a normal length of life, although they are usually completely or par- 
tially sterile. This new allele of W has been designated as W’ by The 
International Committee on Mouse Genetics Nomenclature (Jour. Hered., 
in print). The v stands for viable and since it helps to express the nature of 
the mutant and leaves free for future use the less descriptive symbols W’ 
and W?, it will be used hereafter in this paper. Griineberg (1939) conducted 
haematological studies on different alleles of W, including the new muta- 
tion W’ (W’ of Little and Cloudman) which he designates W?. He found 
that an inherited anemia existed and that the erythrocytes were reduced to 
about two-thirds of the normal values, while the white blood cells were re- 
duced to a lesser extent. He established the normal values by doing blood 
counts on the normal litter mates of the W’W’ animals. 

The following is a report of the results of the microscopic examination of 
some of the organs of 3 male and 3 female mice which were homozygous for 




















VoL. 27, 1941 GENETICS: FEKETE, LITTLE AND CLOUDMAN 115 


the W’ gene. None of these animals had any litters although they were 
mated with their normal sibs. Emphasis, therefore, was placed on studies 
of the sex glands. The ovaries were serial sectioned and several sections 
were cut of the testes. Sections of the accessory male sex glands, and some 
of the endocrine glands, were also examined. Blood counts were taken 
on four of the animals. 

Among the three males No. 7309 was 6 months, No. 5901, 21 months 
and No. 2655, 22 months old when killed. The testes of all of them were in 
normal position but were abnormally small. Microscopically the tunica 
albuginea of all of the glands was unusually thick. The seminiferous 
tubules of the two older animals were lined by Sertolian syncytium and 
seminiferous cells were not present. In the tubules of both testes of the 6 
months old mouse, a few short areas showed spermatogenesis. Groups of 
spermatids and a few spermatozoa were present. In the glands of all three 
animals the interstitial cells were present in larger numbers than normal. 
The efferent ducts and the ductus epididymidis in the caput and corpus 
epididymidis were normal in structure, but contained no spermatozoa. 
The ductus epididymidis in the cauda was distended and contained a 
homogeneous secretion. The ductus deferens was normal. The vesicular 
glands were small, and contained a scant amount of secretion. The coagu- 
lating glands, the ventral and dorsal prostates, and the bulbo urethral 
glands did not show any remarkable deviation from the normal. 

The thymus of 075901 was large and contained several unusually large 
cysts, partially lined by ciliated columnar epithelium. (Such cysts, al- 
though smaller in size, are frequently present in mice.) In the other two 
males this gland was normal in size and structure. The thyroids, para- 
thyroids and adrenals were normal. The pituitaries of all three animals 
were slightly enlarged. Microscopically the anterior lobes showed hy- 
peremia. Although cell counts were not made, the eosinophils (alpha 
cells) seemed to predominate and enlarged vacuolated basophils were not 
present. 

Among the three females, No. 6772 was 11 months; the other two, Nos. 
5588 and 5589, were 22 months old when killed. None of the ovaries of the 
females contained any normal ova. The ovaries of the two older animals 
were composed of dense stroma, granulosa cells and hypertrophied lutein- 
like cells. Many glandular areas lined by simple cuboidal (follicular) 
epithelium were present. In the left ovary of mouse No. 5588 some of 
these areas were cystic. The germinal epithelium showed numerous in- 
vaginations which gave the glands an irregular outline. The left ovary of 
mouse No. 6772 was unusua'ly large and contained a cyst which was filled 
by an organized blood clot. A nodular area composed of densely arranged 
granulosa cells showing many mitotic figures as a sign of rapid growth was 
also present. -In all three animals the epithelial cells of the fimbria and 
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ampulla of each oviduct showed secretory activity and many of the cells 
contained clear vacuoles. Grossly the uterine horns of all three animals 
appeared distended. Microscopically the mucosa of the endometrium 
showed a marked decrease in cellularity and the presence of a large amount 
of hyaline intercellular substance. The endometrial glands were consider- 
ably enlarged and some of them were cystic. 

The adrenal glands of all three females showed some changes in the zona 
glomerulosa which were more pronounced in the older than in the younger 
females. Small groups of cells of this zone were hypertrophied; other small 
densely arranged cells projected in toward the cells of the zona fasciculata 
causing some lobulation, and disturbing the usual orderly arrangement of 
the cells of this zone. The thymus, thyroid and parathyroid glands did 
not show any abnormalities. The pituitaries were slightly enlarged and 
showed hyperemia. 

The skin with the mammary glands in normal position was carefully re- 
moved from each of the three females. Observation with the aid of a dis- 
secting microscope showed that the mammary glands consisted of fairly 
well-developed duct system possessing few side branches. 

Blood counts were carried out on four animals. (The authors are in- 
debted to Dr. L. W. Law and Dr. W. E. Heston for making the blood 
counts.) Table | gives the results. 


TABLE 1 
ToTAL ERYTHROCYTE AND LEUKOCYTE COUNT AND DIFFERENTIAL LEUKOCYTE COUNT 
or W°W® Mice 
DIFFERENTIAL LEUKOCYTES IN % 
NO. OF AGE IN ERYTHROCYTES LEUKOCYTES LYMPHO- MONO- NEUTRO- EOSINO- 
MOUSE MONTHS PER MM.? PER MM.* CYTES CYTES PHILS PHILS 
o'5901 21 4,848,000 3940 72 12 14 2 
9 6772 11 3,264,000 8550 59 12 29 0 
9 5588 22 3,200,000 3200 70 17 13 0 
9 5589 22 3,264,000 Clotted 55 13 32 0 


Average 3,644,000 


The average total erythrocyte count is lower in our animals than the 
figure given by Griineberg. As he correctly states, the blood values vary 
strikingly with age; therefore, comparison of his age group (33-88 days 
with ours (11-22 months) is not advisable. To ascertain the approximate 
differences existing in normal younger and older animals, blood counts 
were taken on C57 black mice. The results show that a reduction of eryth- 
rocyte count can be expected in the older age groups. This reduction is 
comparable in the normal and in the abnormal (W’W’) animals. The re- 
sults are summarized in table 2 and are self-explanatory. 
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TABLE 2 
AVERAGE ERYTHROCYTE Count IN NoRMAL AND W°W’ Mice 

AVERAGE NUMBER OF 

STRAIN OF MICE AGE ERYTHROCYTE COUNT INDIVIDUALS 
C57 black 77-86 days 9,038,000 13 
C57 black 22-23 months 7,636,000 5 
Griineberg normals 33-88 days 9,030,000 8 
Griineberg W°W" 33-88 days 5,937,857 7 
Little W°W’" 11-21 months 3,644,000 4 


Griineberg stated that a reduction of erythrocytes to about one-half of 
the normal value characterizes the homozygous W” animals. A similar 
reduction was found to be present in the W’W” animals reported on in this 


paper. 


Little, C. C., Am. Nat., 49, 727-740 (1915). 

de Aberle, S. B., Am. Jour. Anat., 40, 219-247 (1927). 

Little, C. C., and A. M. Cloudman, Proc. Nat. Acad. Sct., 23, 5385-537 (1937). 
Griineberg, Hans, Genetics, 24, 777-810 (1939). 


OPTICALLY RECIPROCAL GRATINGS AND THEIR 
APPLICATION TO SYNTHESES OF FOURIER SERIES 


By M. J. BUERGER 
DEPARTMENT OF GEOLOGY, MASSACHUSETTS INSTITUTE OF TECHNOLOGY 
Communicated January 2, 1941 


1. Introduction.—Bragg' has recently called attention to the possibility 
of using the diffraction pattern from an appropriately constructed plane 
grating, to give a synthesis of the Fourier series commonly used to display 
the electron density distribution in a crystal. This use of diffraction is 
extremely limited due to the experimental lack of knowledge of the phases 
to assign to the individual holes in the grating, and the difficulty of ex- 
perimentally controlling them when known. The writer? has pointed out 
that an extremely important field of usefulness for the diffraction synthesis 
of Fourier series exists in the syntheses of Patterson and Harker patterns, 
in both of which cases the phases are known. The application is particu- 
larly direct in case of the simple Patterson two-dimensional synthesis, in 
which case all the phases are the same and no special effort is necessary for 
their control. 

Although the diffraction from simple plane gratings is well known, a dis- 
cussion of gratings and their diffraction from a reciprocal lattice point of 
view has not been given. 
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2. Lattice Reciprocity of Plane Gratings.—An important characteristic 
of the diffraction from a plane grating is that the image of this diffraction is 
very approximately a lattice reciprocal to the grating producing it. This 
may be easily proven in the following way: 

Let the original grating be designated G*. (This starred designation is 
used because an important application of this diffraction is to produce a 
crystal lattice, G, from x-ray data in the form of its reciprocal lattice, G*.) 
In figure 1, let a beam of parallelized monochromatic light of wave-length, 
\, strike the grating at normal incidence. Consider, first, a single row of 
the grating (a horizontal row in G*, Fig. 1). All holes in the row scatter 
in phase in a direction contained in a plane at right angles to the row. Now 
consider the conditions under which the several parallel rows will reénforce 
one another’s diffraction within this range of freedom. In figure 2, it is 
evident that this occurs when the rows scatter in directions, 6, to the inci- 
dent beam, such that 

sin 6 = ° (1) 
where d* is the spacing of the rows and m is an integer. If these reénforced 
wavelets are focused by a lens, an image of bright spots occurs at the focal 
plane of the lens. The spots appear in the plane normal to the lattice 
rows in G*, and have positions at distances, ¢, from the center, given by 


t = D tan 6, (2) 


where D is the focal length of the lens. Combining (1) and (2) gives 
r 
t = D tan sin aa (3) 


In actual practice, the spacing of the rows, d*, is easily made of the order 
of a few millimeters, so that this is very large compared with \. Conse- 
quently 6 is very small and a very close approximation is 


sin 6 = tan 6. (4) 


Making this approximation, (3) becomes 
1 a 
t= (Dy)n a* (5) 


The terms in parenthesis are constant for the experiment, hence the several 
t’s are integral multiples of the reciprocal of the row spacing, d*. This 
means that the diffraction due to a set of rows in G* is a lattice row in G 
normal to the original row. The translation interval, ¢, of the new row in 
G is the reciprocal of the spacing of the rows in G*. The image row in G is 
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FIGURE 3 
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thus the reciprocal lattice equivalent of the spacing of the rows in G*. 
Since this same effect occurs for all rows in G* and follows the same law 
as (5) with the same experimental constant, (Dd), the complete diffraction 
of the lattice G* is its reciprocal lattice, G. 

3. Pattern Reciprocity of Plane Gratings.—Reciprocal gratings are re- 
ciprocal in an even wider sense than that of dimensional reciprocity, for the 
pattern of one grating is the pattern reciprocal of that of the other. Figure 
3 shows a section of the gratings of figure 1 taken along the direct ray 
line and normal to any desired row inG*. Let us examine the contributions 
of the several rows of the set in some direction, making an angle, 6, with 
the direct beam, such that the phases of the waves they scatter do not 
necessarily reénforce one another. In this dirction, the zero row scatters 
a wave whose phase may be taken as zero. The first hole away scatters a 
wave in this direction whose phase differs from that of the zero hole by the 
phase angle, y; this is equivalent to a fraction of a cycle equal to ¢/2m and 
to a linear measure of wave amounting to (¢/2r)\. Figure 3 indicates 
that the second row from the origin scatters a wave whose phase differs 
with that of the origin wave by 29, the third hole by 39, etc. Further 
generalizing the conditions of the last section, let the several rows of the 
set have different transmissions: let the transmission of the central row 
be Fo, let the transmission of the first hole be F;, etc. These features of 
the situation may be tabulated: 


PHASE ANGLE OF WAVE FRACTION OF CYCLE 
ROW NUMBER ROW TRANSMISSION SCATTERED BY ROW EQUIVALENT TO PHASE ANGLE 

0 Fo 0 0 
7) 

1 P, ¢ — 
2r 

7) 

2 F; 2¢ 2:— 
; |: : ; lr 
: : : . - 
q Fy qe t= 
2x 


The total contribution of all rows to the wave scattered in the direction, 
5, is the summation of the Fourier series 


S; oi Fe’? + Fee") + Fee? @e/2=)) + Sheen ae | ean (6) 


> F ert i(ale/2)) 
q : 


This summation is performed by the diffraction, and the lens focuses the 
resulting wavelets originating from the rows in G* at a point in G at a dis- 
tance, 7, from the direct beam record, and lying along a line at right angles 
to the rows in G*. This distance, r, can be computed according to (2) as 


r = D tan 6, (7) 
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and the direction angle, 4, is given by a relation similar to (1), except that 
the integer, , is replaced by a cycle fraction, ¢/27: 


7) 
——) 
sin 6 = =. (8) 
Combining (7) and (8) gives 
xh 
r = D tan sin- _* (9) 


Making approximation (4), (10) becomes 


2a 
= D —., 
r 7 (10) 
Dividing (10) by (5), and letting » = 1, gives 
r ¢g 
” ise Qa (11) 


This means that the cycle fraction, ¢/27, may be replaced by a fraction, 
r/t, which represents the distance, 7, as a fraction of the full translation, 
reciprocal to the row spacing in question. Making this replacement in (6), 
the summation becomes one pertaining to the point at a distance, 7, in G, 
rather than a summation for the direction, 5, from the normal to G*: 


S, = bated (12) 


The diffraction due to the set of rows selected in G* occurs not only at 
the single point having a central distance, 7, in G, but at a continuum of 
points embracing all values of 7, and forming a radial streak through the 
center of G. The full diffraction of the grating G* comes not only from this 
particular row, but from all conceivable rows, no matter how complex. 
In this fashion, every radial streak in G is the pattern reciprocal of some 
conceivable set of rows in G*. 

4. Pattern Synthesis in Terms of Lattice Coérdinates.—It is desirable to 
recast the results of the last section in terms of coérdinates referred to lat- 
tice axes chosen in G* and their reciprocals in G. If axes b* and a* are 
chosen in G*, the diffraction due to sets of rows parallel with them pro- 
duces radial streaks in G parallel to a and ), respectively. The summations 
for the particular point at a distance, x, along a and for the particular 
point at a distance, y, along b, are, according to (12), 
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Ss ie THe eee) (13) 
wind S, = TRereo (14) 


We now require to know how each hole in G*, as indexed on the axes, 
a*, b*, contributes to the summation at a point having absolute linear 
coérdinates, x, y, with respect to the axes, a, 6, in G. Figure 4 shows 
how a hole in G*, having index hk, contributes to diffraction in some gen- 
eral direction. In this direction, the holes in neighboring positions along 
the a* axis scatter with a difference in phase represented by the fraction 
of a cycle, g,*/27. The Ath hole along the a* axis therefore scatters a 
wave whose phase differs by /(¢,+/27) cycles from that of the origin hole. 
Similarly, the holes in neighboring positions along the b* axis scatter with a 
difference in phase represented by the fraction of a cycle, g»«/27. The 
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FIGURE 4 


kth hole along the b* axis therefore scatters a wave whose phase differs by 
k(gp+/2mr) cycles from that of the origin hole. The total phase difference 
between the hole, /k, and the origin hole is therefore 





a + ee. (15) 
Qa 2Qr 


Pie =h 


From relations similar to (11), 
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Substituting these in (15), the phase of the hole, hk, with respect to that of 
the origin, may be very closely approximated by 


Pip © he he (18) 
a b 


Now this wave scattered by the hole, hk, is focused at the point having 
coérdinates x, yin G. That this is the case can be easily grasped from the 
following reasoning: In (15), let y«/2m successively take the particular 
constant values of 0, y/b, 1. For each of these values, let y,+/2r vary a 
little from zero. In this variation, the focused ray traces out a hyperbola 
on G because the locus of ¢,+/2r = constant is the surface of a cone which 
intersects G in a hyperbola. Since (gj«/2m)\ is exceedingly small com- 
pared with a* in practical gratings, this hyperbola is substantially a straight 
line, a condition correlating with the approximation of (4). These several 
lines are labelled g+/2r = 0, ¢«/2r = y/b and ¢«/2r = 1 in figure 4. 
In the same way, ¢,+/2r = 0, g+/2r = x/a and ¢+/2r = 1 may be 
mapped in figure 4. From this it follows that the wave having the phase 
indicated by (15) is focused at a point in G having relative coérdinates 
x/a, y/b (i.e., at absolute coérdinates, x, y). There is an inconsequential 
approximation involved here, but the approximation in no way affects the 
synthesis, but merely the position where the synthesis is focused. 

From this discussion, it follows that at the point in the grating, G, having 
coordinates x, ‘y with respect to the grating, axes, there occurs the synthesis, 


Si, por DF alla + R(y/d)) (19) 


This is in exactly the same form as the usual expression for the Fourier 
synthesis of the electron density at the point, x, y ,in the two-dimensional 
projection of a crystal cell. The corresponding Patterson synthesis can be 
obtained by making the holes in grating, G*, proportional to |F|* rather 
than F, giving the summation 


Mig: YD] Fal 2 oe’ + R(y/d)} (20) 
which may also be written 


Any * LD Frel? cose” (4@/a) +k/0)} (21) 


5. Control of Phases.—In electron density syntheses, the transmission, 
Fiz, is, in general, a complex quantity usually of unknown phase. Two 
methods have been suggested” ? for the control of this phase for the cases 
where the phase angle is known to be 0 or 7. Any phase may be impressed 
on a term by placing in the path of each grating hole a thin sheet of plane- 
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parallel transparent material whose angle with respect to the transmitted 
beam may be varied, thus varying its effective optical thickness. Pieces 
cut from a single, perfect, mica cleavage are suggested for this purpose. 

Another device useful in the control of phases is the utilization of mul- 
tiple parallel gratings. Each gives a separate diffraction pattern, and 
these patterns may be optically combined, for example, with the aid of 
mirrors. The phase of an entire grating may be controlled, as in the last 
paragraph, with the aid of a single sheet. The use of two gratings having 
different phases may be applied to certain Harker syntheses, in which 
positive and negative terms must be summed at the same point of the 
pattern. The use of two (or more) gratings having different phases may 
also be applied to electron density summations, by permitting transmission 
of a given index, hk, to be contributed only from a single grating of the 
several used. For example, for centrosymmetrical projections, one grat- 
ing may be miiie to have a positive phase, the other negative. A hole, hk, 
is uncovered on only one of the two gratings. This method may be ex- 
tended to non-centrosymmetrical projections, by using several gratings; 
one grating transmits for the symmetrically equivalent /,k;’s, another for 
heke’s, etc. In this case, each grating has a general phase. 

1 Bragg, W. L., ‘““A New Type of ‘X-Ray Microscope,’ ’’ Nature, 143, 678 (1939). 


2 Buerger, M. J., ‘“The Photography of Interatomic Distances Vectors and of Crystal 
Patterns,”’ Proc. Nat. Acad. Sci., 25, 383-388 (1939). 
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SEROLOGICAL REACTIONS WITH SIMPLE SUBSTANCES 
CONTAINING TWO OR MORE HAPTENIC GROUPS 


By Linus PAuLinc, DAN H. CAMPBELL,* AND DAviID PRESSMAN 


GATES AND CRELLIN LABORATORIES OF CHEMISTRY, CALIFORNIA INSTITUTE OF 
TECHNOLOGY 


Communicated January 13, 1941 


It has been inferred from inhibition studies that a simple hapten can 
combine with antibodies homologous to a conjugated protein involving 
that hapten! to give soluble compounds. In most cases haptens do not 
give the precipitin reaction with antiserum nor produce anaphylactic shock 
in animals sensitized with corresponding conjugated proteins. It was 
pointed out recently® that the framework theory of antibody-antigen pre- 
cipitates suggests strongly that simple substances containing two or more 
haptenic groups should, under favorable circumstances, show all of the 
properties of more complex antigens; the thorough experimental test of this 
should provide evidence in support of or against the framework theory. 


Landsteiner and van der Scheer* observed the occurrence of the precipitin 
reaction and anaphylaxis with simple substances formed by coupling two 
anilic acid molecules with resorcinol or tyrosine. On the other hand, 
Hooker and Boyd‘ have reported that they observed no forma- 
tion of precipitate with several different compounds containing two hap- 
tenic groups. In order to obtain further evidence, we have prepared 
several substances with two or more haptenic groups per molecule and have 
tested them with suitable serum and with sensitized guinea pigs. We ob- 
served precipitate formation with each of the substances, and also observed 
anaphylaxis and desensitization with one of them, thus verifying the 
existence of the phenomena reported by Landsteiner and van der Scheer. 

The serum used in our experiments was obtained from rabbits which had 
been given a course of six injections (alternately intraperitoneal and in- 
travenous) of azoprotein per week for four weeks, the rabbits being bled 
nine days after the last injection. The azoprotein used was made by 
coupling sheep serum protein with diazotized arsanilic acid. Four azo- 
protein preparations, containing 0.68%, 1.9%, 2.3% and 3.3% arsenic 
(relative to total protein as determined from Kjeldahl nitrogen), respec- 
tively, were used. The titers of the serum of fourteen rabbits given one or 
another of these preparations were nearly the same, about 0.4 mg. anti- 
body per ml., as found by Kjeldahl nitrogen analysis of the precipitate 
given with azoproteins made from egg albumin and arsanilic acid. 

The haptens prepared and tested and the results obtained are shown in 
table 1. 
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VII 
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TABLE 1 


AMOUNT OF PRECIPITATE 


None 


Small 


Medium 


Medium 


Small 


Medium 


Large 


R’ = —NN¢ NNC > AsOsH 


The haptens I to VI were made by reaction of phenol, resorcinol, phloro- 
glucinol, or p,p’-dihydroxybiphenyl with diazotized arsanilic acid, and VII 
by coupling with p-nitrosophenylarsonic acid the triamine obtained by the 
hydrolysis of the product of the reaction of phloroglucinol with diazotized p- 
The haptens were dissolved in alkali, brought to pH 
between 7 and 9, diluted and mixed with serum (0.5 ml. hapten solution to 
At suitable concentrations of each of the haptens II to 
VII, between 10~‘ and 10~, a noticeable precipitate could be seen in a few 
minutes, becoming heavier on standing overnight in the refrigerator. No 
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precipitate was obtained in control tests made with normal serum, and no 
precipitate was given by hapten I with antiserum. The optimum hapten 
concentration was in each case about 10-5 (0.01 mg. hapten per ml.). 

Well-defined precipitates were obtained with substances II to VI, with 
no obvious dependence on variation of the number of haptenic groups 
from two to four. Substance VII gave a better reaction than the others; 
this we attribute to the smaller steric repulsion of antibody molecules 
attached by the long azobenzeneazophenylarsonate groups R’ than that 
of those attached by the shorter azophenylarsonate groups R. 

The amount of antibody precipitated by substance VII under optimum 
conditions, about 0.4 mg. per ml. of serum, was nearly as great as that 
precipitated by azoprotein. We did not analyze the precipitates for hap- 
ten; with the assumption that all the hapten present is precipitated in the 
antibody-excess zone, our antibody analyses give for the antibody/antigen 
molecular ratio values between 0.2 and 0.9. These are smaller than the 
values found for protein antigens. 

Guinea pigs actively sensitized by injection of azoprotein or passively 
sensitized by injection of rabbit anti-arsanilate serum were found to suffer 
shock on intravenous injection of '/2 ml. of solution containing 1 to 4 
mg. of hapten VII, and to be desensitized with respect to azoprotein. Pre- 
liminary efforts to produce active antisera in rabbits and to sensitize guinea 
pigs by injection of hapten solutions gave negative results. 

Our investigations confirm those of Landsteiner and van der Scheer. 
We attribute the antigenic properties of the substances studied to the pres- 
ence in their molecules of two or more haptenic groups, and we consider 
that the phenomena provide strong support for the framework theory of 
serological precipitates, which provides a simple explanation of them. 

Landsteiner® has suggested as a possible alternative explanation that 
“the ready precipitability of these dyes is dependent upon peculiarities in 
constitution which, like those of fatty groups, diminish solubility in water 
and favor the formation of colloidal solutions.” Landsteiner and van der 
Scheer reported somewhat better precipitations with hapten solutions 
which had aged by standing than with freshly prepared solutions. We, 
however, have observed no difference in behavior of old and fresh solutions. 
Our haptens dissolved easily at pH 9 to give clear stable solutions. At pH 7 
the solutions precipitated on standing overnight, but mixtures with normal 
serum remainéd clear. 

The failure of Hooker and Boyd to obtain precipitates with simple sub- 
stances containing two haptenic groups may be attributed in part to the 
small size of the molecules of some of the substances studied by them, such 
as biphenyl-p,p’-dicarboxylic acid, 


HOOC< >< coon 
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and p,p’-dicarboxylazobenzene, 


HOOCK —_NNC SCOOH. 


Steric repulsion between two antibody molecules attached to such small 
molecules would be much stronger than for the molecules used by Land- 
steiner and van der Scheer and by us. This explanation does not apply to 


OH OH 


HOOC’ oy NNC COOH, 


HOS A S0;H 


for which they also report negative results. 

Summary.—We have verified the observations of Landsteiner and van 
der Scheer that certain simple substances (azodyes) give precipitates with 
antisera and produce anaphylactic shock in sensitized animals. The sub- 
stances found effective by these authors and by us contain two or more 
haptenic groups. The difference in properties of these substances and 
those containing only one haptenic group, which form only soluble com- 
pounds with antibodies, is explained in a simple way by the framework 
theory of serological precipitates, and the observations provide strong 
support for this theory. 


* Fellow of the Rockefeller Foundation, 1940. 

1 Landsteiner, K., The Specificity of Serological Reactions, Charles C. Thomas, Balti- 
more (1936), Chap. V; Biochem. Z., 93, 117 (1919); 104, 280 (1920); Marrack, J., and 
Smith, F. C., Nature, 128, 1077 (1931); British J. Exp. Path., 13, 394 (1932); Hauro- 
witz, F., and Breinl, F., Z. phystol. Chem., 214, 111 (1933). 

2 Pauling, L., Jour. Am. Chem. Soc., 62, 2643 (1940). 

? Landsteiner, K., and van der Scheer, J., Proc. Soc. Exp. Biol. and Med., 29, 747 
(1932); Jour. Exp. Med., 56, 399 (1932); 57, 633 (1933); 67, 79 (1938). 

* Hooker, S. B., and Boyd, W. C., Proc. Third Internat. Cong. Micrcbiol., p. 814 
(1940); personal communication. 

5 We acknowledge the assistance of Mr. Miyoshi Ikawa and Mr. Carol Ikeda in the 
preparation of these compounds. 

§ Landsteiner, K., The Specificity of Serological Reactions, p. 120. 
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PROGRESSIVE IMMUNIZATION 
By Epwin B. WILSON AND JANE WORCESTER 
HARVARD SCHOOL OF PuBLIC HEALTH 


Communicated January 3, 1941 


In a recent note! we discussed the immunization of a population through 
contact of its individuals with a virus on the assumption that at the first 
contact the disease was acquired with a resulting lasting immunity. Such 
an all-or-none hypothesis with respect to the process of individual immun- 
ization may be sufficiently realistic for a disease like measles. Conceiv- 
ably, however, a population might gain in immunity to some diseases 
through its individuals gaining progressively in immunity by successive 
contacts with the pathogenic agent—indeed as an extreme case it might be 
imagined that complete immunization was obtainable in that way without 
the development of cases of the disease.” 

Suppose that if a person who has had no previous contact with the 
pathogen comes in contact with it, he has a chance /; of contracting the 
disease and acquiring a complete immunity and a chance g, = 1 — p; of not 
acquiring the disease but of gaining some immunity. When the person 
who has not had the disease comes a second time in contact with the 
pathogen he will have a chance p2 of acquiring the disease with ensuing 
complete immunity and a chance q2 of escaping but with the possibility of 
gaining somewhat further in immunity. The fractions of the persons just 
after the second contact will then be three: 


pi who acquired complete immunity on first contact. 
dip, who acquired complete immunity on second contact. 
9:92 who remain not completely immune. 


Continuing the argument one finds four fractions after the third contact, 
namely, ~1, giPe, 91923 who have had the disease and g:g2g3 who have been 
thrice in contact without acquiring it but presumably have gained in im- 
munity. 

By the law of small numbers, if N contacts per capita have been spread 
over a population the fraction of persons who have had just k is e~™ N*/k!. 
If there be dN further contacts per capita, the fraction of persons who have 
their (k + 1)s¢ contact is dN multiplied by the previous fraction, and of 
these the fraction igo . . . d¢ P41 Will contract the disease and gain com- 
plete immunity. If C(x)dx be the case rate during the time dx when dN 
contacts are distributed, we have*® 


k= @ 


a Nt 
C(x)dx = e "2 auae ~ + + UePeti ey = o(N)dN. 
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The average number N of contacts which have been experienced is itself a 
function of age x, but we need not make any special assumption with 
reference to it. The total number of cases per capita up to age x is 


Y(N) = ff cons = (1 — Q)(1 — e~*) 


. r N? N3 
—e ‘| @-ON + (qge—O) + (Gigegs — Qa + - - | 


where (Q stands for the indefinitely continued product? of the q’s. 
In no actual case do we know the values of p;, but an assumption will 
be made merely in way of illustration.’ Suppose 


/ / k 
Pr _ 1/9, pe = 45 ae » Pr sacs V/s 


OSs eee ye 
Q= 5°75: oe + converges to 2 = 0.2888. 


Thus after an infinite number of contacts (V = ©) the total fraction of the 
population which has acquired the disease is 1 — Q = 0.7112 and the frac- 
tion Q = 0.2888 have acquired immunity progressively by repeated contact 
with the pathogen without ever contracting the disease. The values of 
¢(N), ¥(N) and the increment Ay(JV) for successive integral values of N 
from 0 to 7 are given in table 1. If the lapse of time (or age) be measured 
by the increase of the average number of contacts N, the column ¢(V) may 
be called the instantaneous case rate (per capita) by virtue of the relation 
C(x) = o(N)dN/dx, the value of ¥(V) may be called the total accumulated 
case rate, and the increments Ay(JV) are the case rates for the intervals from 
N to N + 1 measured by average contacts. 

It may be recalled that on the all-or-none hypothesis, the corresponding 
formulas! are 


o(N) = e~%, WN) = 1 — e7%, Av(N) = eX(1 — 1/e). 


Here log ¢(N) and log Ay(N) plot as linear functions of N and ¥(N) is a 
saturation curve which gives the fraction J(V) of the population which has 
acquired immunity through having the disease. If we be dealing with the 
all-or-none hypothesis but with the complication, which is surely present in 
practice, that only a fraction of the cases are reported, whether this be 
current reporting of communicable diseases or results found by surveys 
which ask whether persons have had the disease,* and if we assume for 
simplicity that the under-reporting is independent of age (or of NV) we 
should have to replace these formulas by 


g(N) = fe~%, WN) = f(1 — e~%), AY(N) = fe~*(1 — 1/e) 
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in terms of the supposed true condition in the population; such under- 
reporting could not fail to be present in respect to sub-clinical cases even if 
there were none with respect to the clinical cases. The only effect on our 
graphs would be that log,g would be — N + log f in place of log, e = —N, 
and similarly for log Ay; the lines would still be straight, but below their 
true positions when based on reported cases instead of on all cases. 

If the values of log ¢ or of log Ay in table 1, as computed from our 
special hypothesis relative to progressive immunization of individuals, be 
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Graphs of logio ¢(N), solid dots, and logio Ay(V), crosses, plotted against VN. The 
upper graph I is for table 1 and the lower graph II is for table 2. 


graphed, it will be observed that these lines are not far from straight (Fig. 1, 
I). Indeed with the inaccuracy of observed material what it is at present, 
even at best, the departures of the observations from any assumed law un- 
doubtedly would be greater. Thus we may assert that if the population 
were immunized progressively as assumed, the observations could fairly be 
interpreted as due to the all-or-none hypothesis with fewer contacts and 
with under-reporting. Indeed the formulas 


y(N) = 0.4635e7°""-¥(N) = 0.7193(1 — e~**# 9), 
Ay(N) = 0.7193e7°* Nai — e7 omy 
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fit the data in table 1 very well. If N’ = 0.6444 N, so that the contacts 
come only 0.6444 times as fast, and if the disease be considered only 71.93% 
reported, the results of the all-or-none hypothesis would practically dupli- 
cate’ those in table 1. 

Suppose in the second place that 


pi = 0.04, pe = 0.03, ps = 0.02, ps = 0.01, p; = 0,1 > 4 


so that there is never much chance of a case arising from contact with the 
pathogen, and no chance at all after four contacts* (because of the progres- 
sive acquirement of complete immunity). In table 2 are given the values 
of gy, y, Aw and in figure 1, II, are graphed log g and log Ay. The graphs 
are again not far from straight. The table could be explained as arising on 
an all-or-none hypothesis with NV’ = 0.5165 N, i.e., with the contacts com- 
ing only about half as frequently, and with a reporting of only about ten 
per cent, whether because of failure to report clinical cases or because of the 
large number of sub-clinical cases which gave immunity.° 


TABLE 1 TABLE 2 

N ¢ v A ¢ ¥ Ay 
0 0.5000 0.000 0.0400 0.0000 

0.353 0.0346 
1 0.2400 0.353 0.0293 0.0346 

0.172 0.0245 
2 0.1190 0.525 0.0199 0.0590 

0.086 0.0160 
3 0.0609 0.612 0.0125 0.0750 

0.045 0.0098 
4 0.0321 0.657 0.0074 0.0848 

0.024 0.0056 
5 0.0174 0.681 0.0041 0.0904 

0.013 0.0030 
6 0.0097 0.694 0.0022 0.0934 

0.007 0.0016 
7 0.0055 0.701 0.0011 0.0950 


For the case of this second example we give in table 3 against the average 
number of contacts, for the total population the number of persons per 
thousand who have had 0, 1, 2, 3 or 4 or more contacts and (in italics) the 
number who have had the disease. By our assumptions, those who have 
had four or more contacts are immune whereas those who have fewer con- 
tacts and have not had the disease are not yet completely immune. 

These two specific hypothetical illustrations in which we have seen that 
progressive immunization simulates that of the all-or-none hypothesis with 
partial reporting are offered as a comment on our remark’ that it is not 
necessary to infer from the similarity of age distribution of two diseases of 
childhood that the immunizing process relative to the individual is similar ; 
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TABLE 3 


NUMBER OF PERSONS PER THOUSAND WHO HAVE Hap A GIVEN NUMBER OF CONTACTS 
AND (IN ITALICS) NUMBER WHO HaAvE Hap THE DISEASE 


AVG. ACTUAL NUMBER OF CONTACTS TOTAL OTHER 
N 0 1 2 3 4+ CASES IMMUNES 

1 0 15 13 5 2 35 17 
368 368 184 61 19 per thousand 

9 0 11 19 16 14 59 129 
135 271 271 180 143 per thousand 

3 0 6 15 20 34 75 319 
: 50 149 224 224 353 per thousand 

4 0 g 10 17 55 85 512 
18 73 147 195 567 per thousand 

5 0 1 6 12 71 90 664 
Y 6 34 84 140 735 per thousand 

6 0 1 3 8 82 93 767 
2 15 45 89 849 per thousand 

~ 0 0 2 § 8&9 95 830 
1 6 22 52 918 per thousand 


for one disease it may be all-or-none, for another it may be progressive with 
essentially similar curves of age distribution. Considerations other than 
these have to be brought to bear to make the decision—at least until the 
available material is decidedly better. As the all-or-none hypothesis is the 
simplest, one may appeal to the principle of Occam’s razor to hold fast to it 
except as definite immunological observations require us to depart from it.” 


1 These PROCEEDINGS, 27, 7-13 (1941). The assumption was made, and is continued 
in the present article, that the mortality due to the disease is negligible. 

2 For the purpose of this discussion a ‘‘case’”’ corresponds to the completion of the 
acquisition of immunity all at once whether the case be clinical or sub-clinical and thus 
contrasts with a gain in immunity which remains only partial. 

3 Calculations will be made in terms of the average number of contacts and not in 
terms of age because of the added complication which arises from the uneven distribution 
of contacts by age discussed in our previous note.! The use of the law of small numbers 
implies a population homogeneous with respect to exposure except for the operation of 
chance. 

4 The product Q = (1 — p:)(1 — po)... (1 — p,) ... will converge to a value between 
0 and 1 if the series p: + po... + py +... converges, but will converge to 0 if the 
series diverges. Thus if we had p; = !/2, po = 1/3,..., Dy = 1/k, . . . we should have 
Q = 0 and everybody who lived long enough would contract the disease despite any 
immunity previously acquired by contact with the pathogen, the evidence of progressive 
immunity being merely the progressively lessened likelihood of getting the disease on 
successive contacts. 

5 The aim of this theory of progressive immunization with its illustration by specific 
examples is like that in Zinsser and Wilson, ‘‘Bacterial Dissociation and a Theory of the 
Rise and Fall of Epidemic Waves,” Jour. Preventive Med., 6, 497-514 (1932), namely, 
to show the consequences of what might happen rather than to maintain that such a 
thing does happen. , 

6 The factors of under-reporting may not be the same in the two cases; for example, 
surveys of persons for a history of measles generally give values of about 90% at ages 15 
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and over, which cannot be more than 10% under full reporting, whereas reports of cases 
to departments of public health are very much more deficient. 
7 If we assume p, = ar*—,a < 1,7 < 1, then 
(rN)* 
3! 





2 
¢g(N) = ae~*1 1 + qrN + ang + 919293 Te ]<eonve" 


since all g’s are less than 1; but all products of g’s are greater than Q, hence ¢(N) > 
aQe~Ne"’. Hence if @ be an appropriate fraction 


e(N) = ae—?-*%[1 — 0(1 — Q)]. 


Now if Q be near 1, i.e., if not many persons contract the disease (finish acquiring their 
immunity all at once), g(N) takes a form interpretable on the all-or-none hypothesis as 
due to contacts coming only 1/(1 — 17) times as frequently with only the fraction a of 
cases recorded as such. 

8 Although four contacts completely immunize by hypothesis we do not consider all 
fourth contacts of persons who have not previously had the disease to be cases, but only 
1% of them; in this there is a verbal inconsistency with the definition laid down in foot- 
note 2. 

® In neither of these illustrative examples is the relationship between log ¢(N) or 
log Ay(N) and N strictly linear. It is possible to choose the successive values ~; in 
such a manner that the relationship is strictly linear, namely, pi+: = apfi/qi with a any 
number less than g:. In this case 

’ ft has —N: * ’ 
V(N’) = i a — e—"’) with N’ = (1 — a)N 
and hence the case rates will be strictly as on the all-or-none hypothesis with contacts 
only 1 — a@ times as frequent but with the cases reported only ~,:(1 — a) of the total 
which arise. For measles, if we assume 90% reporting (as on surveys) we have 
fis(1 — a) = 0.9 and any values of p, and a which satisfy the linear equation 10p, = 
9 — 9a with a < 0.9, will be satisfactory. These are such as 


a 0 0.2 0.5 0.8 0.9 

hi 9 0.72 0.45 0.18 0.090 
bs 0 0.51 0.41 0.18 0.089 
bs 0 0.21 0.35 0.17 0.088 
bs 0 0.05 0.26 0.16 0.087 


and soon. The first column is “all-or-none” with 90% reporting; the others are “‘pro- 
gressive” at various rates but also with 90% reporting. 

If we had to rely on public health statistics where the reporting (though actually vary- 
ing markedly with age) usually averages no more than 50% and if for illustrative pur- 
poses we use the average figure of 50% at all ages we have 


a 0 0.2 0.5 0.8 0.9 

fi .50 0.40 0.25 0.100 0.050 
pr 0 0.133 0.167 0.089 0.047 
Ps 0 0.031 0.100 0.078 0.045 
Ps 0 0.006 0.056 0.068 0.042 


and soon. The results for a = 0.9 imply frequent contacts (ten times as many as on 
the all-or-none hypothesis) but low and very slowly diminishing attack rates upon con- 
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tact, whereas those for a = 0.02 imply only 25% more frequent contacts than the all-or- 
none hypothesis but with rather high though rapidly diminishing attack rates. 

10 Schuman and Doull, Amer. Jour. Pub. Health, 30, Supplement to March, 1940, p. 21, 
state: ‘From these estimates of carrier prevalence and from the average annual in- 
crement in Schick-negatives, an estimate may be made of the number of carrier infec- 
tions necessary, on the average, to produce sufficient antitoxin to give a negative Shick 
test.” This seems to imply belief on their part that immunization as indicated by a 
negative Schick test may be acquired progressively; it is possible that a degree of im- 
munity sufficient to prevent the disease arising on further contact is acquired before the 
Shick test turns negative, and if this degree of immunity were regularly acquired on the 
first carrier infection, we should classify the immunity as arising from a sub-clinical case 
under the definitions of footnote 2 and the phenomenon as to be considered of the all-or- 
none type. 


ON DIRICHLET’S DIVISOR PROBLEM 
By AUREL WINTNER 
DEPARTMENT OF MATHEMATICS, THE JOHNS Hopkins UNIVERSITY 
Communicated January 8, 1941 


Let d(m) denote the number of divisors of the positive integer m, and let 
D(x) = d(1) + d(2) +... + d([x]). 


It is an elementary fact that this summatory function can approximately 
be represented by 


F(x) = x(log x + 2C — 1) 


as x — o, where C is the Euler-Mascheroni constant. The classical 
“divisor problem of Dirichlet’’ deals with the remainder term, D(x) — F(x), 
of this elementary approximation. Dirichlet proved that 


— 


D(x) = f(x) + Ox’), | | # > ©, 


holds for c = !/3; as indicated in his correspondence with Kronecker, he 
has known a sharper estimate also.1 Today? it is known that the O-term is 
valid for a certain c < 1/3, and that it cannot be true for c = 1/, (although 
it might hold for every c > '/,). The present note leaves aside this difficult 
problem of estimation, and concerns itself with the fascinating fact that 
the remainder term, D(x) — F(x), admits of an anharmonic analysis (on a 
proper scale). 

Let Z(s) be an abbreviation for the difference Yi(s) — Ay(s), where VY, 
H, denote the usual transcendents of the theory of cylinder functions. Ac- 
cording to Voronoi,’ the convergent explicit representation 
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D(x) = F(x) + 1/4 + x” - n~""d(n)Z(4en'*x'/*) 
n=1 


_ holds for every x, if D(x) is defined as the limit of !/2{D(x + h) + D(x—h)} 
‘ as h tends to zero. On applying’ to*each of the terms Z = Y, — A, 
the asymptotic formulae of the cylinder functions Y;, H;, one easily obtains 
the convergent expansion 


a x/* © d(n) t/a t/ al 
D(x) = F(x) + V/a+ ah :¥ “7 COS (4an’x’* — 2/4) + O(x”); 


(it is even known that the trigonometric series on the right is boundedly 
convergent on every finite x-range). 

It seems to be natural to replace the independent variable x by ¢t = x’ 
and to ask whether or not the resulting trigonometric series in t is the Fourier 
series of the function to which it belongs in a formal fashion; that is, of the 
function 


/s 
’ 


Q(t) = 1D(#) — F(e)}/t”. 
It turns out the answer is affirmative, in the sense that the reduced remainder 
term, Q(t), of Dirichlet’s divisor problem is almost periodic (B*), 1 <t< ~, 
and has the Fourier series 


ae ‘ 
Q(t) ~ avy a se cos (4arn hy = 1/4). 


This implies, in particular, that Q(¢) has an asymptotic distribution func- 
tion. 

It is clear from Bessel’s identity that, if /} f(t)} is an abbreviation for 
the average 

SR ee 
To FY Sods, 
the italicized theorem is equivalent to the following pair of assertions: 

(i) The amplitude belonging to an arbitrary frequency, that is, the 
Fourier average M{e™Q(t)}, exists for every real number \, and is zero 
or has the value suggested by Voronoi’s formal series according as ) is not 
or is of the form + 41n’/*, where n = 1, 2,.... 

(ii) The average intensity, /{Q(t)*}, exists and is equal to the square 
sum of the amplitudes. 

Cramér* deduced from the formal series a mean-value theorem which, 
in the above notations, may be expressed by saying that 


fed 2 ee 
Tool (4) dt 
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exists and has a value compatible with (ii). Two partial integrations show 
that this result of Cramér is an Abelian consequence of (ii); in fact, 


. ‘ T 
MIQ0} = po, pS Oral 


Thus the converse inference is Tauberian in nature. On the other hand, 
it seems to be hard to find a Tauberian theorem applicable to the present 
case. Consequently, a direct approach is needed. 

Correspondingly, the proof of (i)-(ii), though straightforward, is some- 
what lengthy, and will be published elsewhere. It applies also when 
Dirichlet’s divisor problem, which can be thought of as the problem of lat- 
tice points in hyperbolas, is replaced by the problem of lattice points in 
circles. 


1 Lejeune-Dirichlet, P. G., Werke, 2, 49-66 and 407 (1897). 

2 For references cf., e.g., Enc. der Math. Wiss., II3, 816-820. 

3G. Voronoi, Ann. Ec. Norm. Sup., ser. 3, 21, 207-268 and 459-534 (1904). 
4H. Cramér, Math. Ztschr., 15, 201-210 (1922). 


UNIQUENESS THEOREMS FOR CONFORMAL MAPPING OF 
MULTIPLY CONNECTED DOMAINS 


By M. SHIFFMAN 
DEPARTMENT OF MATHEMATICS, COLLEGE OF THE City OF NEW YORK 
Communicated December 31, 1940 


In a previous issue of these PROCEEDINGS,' it was proved on the basis of 
the Plateau problem that any k-fold connected domain is conformally 
equivalent to one of the following class of domains: The exterior of k non- 
intersecting closed curves which are, respectively, homothetic to k given 
convex analytic curves C;, C2, ..., C,. It was furthermore possible to 
consider these domains normalized, e.g., for k = 3, by varying three of the 
boundaries Ci, C2, C; only by radial transformations relative to a fixed in- 
terior point of each. We shall prove here that a k-fold connected domain 
can be mapped conformally onto a unique one of this class of normalized 
domains. 

It suffices to prove the following: if f(z) is an analytic function mapping 
one domain G of our class of normalized domains conformally on another 
normalized domain G’, then f(z) =z. The procedure is to suppose that 
the analytic function 


g(z) = fe) — 
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is not identically zero (and consequently not identically constant), and to 
obtain a contradiction by counting the zeros and poles of g(z).? 

Let Ci, C2, ..., C, be the boundaries of G, and C,’, C’,..., C,’ of G’. 
Let N, P be the number of zeros‘and poles, respectively, of g(z). It is 
easy to see that P < 2. For, either f(z) takes ~ into a finite point and 
finite point into , in which case P = 2; or f() = ©, for which P = 1 
or 0. Concerning the zeros of g(z) we have 


k 1 7 k 
y-f. TT = >], 
j=l 271 Cj °8 g(2) j=l 7 


where the integration around C; is taken in the clockwise direction, and 
where any zeros of g(z) which occur on C; are circumvented by protrusions 
outside Cj. 

Now, f(z) maps C; into the curve C;’, obtained from C; by a homothetic 
transformation; and as z moves clockwise once around C;, the image point 
f(z) moves clockwise once around C;’._ The complex number g(z) = f(z) — 
z on C; is the vector with initial point z on C; and end-point f(z) on C;’. 
The vectors g(z) form a vector field defined on the curve C;. And if g(z) ¥ 
0 on C;, I; is exactly the order or characteristic of this vector field. We shall 
prove that J; < —1forj = 1, 2,3, and J; < Oforj > 3. 

Consider J;. Since C,’ is obtained from C; by a radial transformation - 
relative to an interior point 0, and since C; is convex,’ it follows that C,’ 
either surrounds Cj, or is interior to Ci, or coincides with C;. In the first 
case, J, = —1. To prove this, contract the curve C,; continuously till it 
approaches a point interior to C,’. In this contraction, each vector re- 
mains attached to its position on C; and on C,’. Since no vector ever be- 
comes zero, the characteristic remains constant. In the limit the charac- 


teristic is —1, so that originally J; = —1. 
A similar result applies if C,’ is interior to C,. Finally, suppose that C,’ 
coincides with C;. If g(z) has no zeros on Ci, expand C;’; i = —1. If 


g(z) has zeros on Cj, the circumvention of a zero will contribute an amount 
—m/2 to the integral J;, where m is the order of the zero. Choose this cir- 
cumvention to be an arc of a circle K of such small radius that the integral J; 
around it differs from —m/2 by less than e, and such that g(z) at each end- 
point of K differs from the tangential direction to C; by an angle less than e. 
On an arc of C, between two successive zeros of g(z), the change in angle 
differs by less than 2e from the change in angle of the tangential rays on Cy. 
Thus, J; differs from —1 — n/2, where n is the sum of the orders of the 
zeros of g(z) on C;, by at most a fixed multiple of e. Hence J; = —1 — 1/2. 

In all cases for the integral J,, therefore, J; < —1. Similarly, 2 < —1 

Consider I; for] > 3. There are three possible situations of C;’ relative 
to C; (besides the possibilities enumerated above for C;, which would yield 
I j <- 1) : 
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1. C;,’ and its interior lie completely outside C;. The characteristic J; 
of the vector field g(z) is 0. 

2. C;’ is tangent to C;. If g(z) has no zero on C;, then C;’ can be dis- 
placed slightly so as not to intersect C;. Hence J; = 0 or —1 according as 
the tangency is external or internal. If g(z) has a zero on C;, a simple dis- 
cussion of the circumvention of this zero shows that J; < —1. 

3. C;’ intersects C; in two points, P and Q.‘ Denote the arcs of C; 
interior and exterior to C;’ by a, B, respectively, and the arcs of C;’ interior 
and exterior to C; by y, 6. Let f(z) map P, Q into the points P’, Q’ of C;,’, 
and the arc ainto a’. Suppose first that g(z) has no zeros on C;, so that P’, 
Q’ are distinct from P, Q, respectively. We shail distinguish three cases: 
either a’ Cd or a’ D 5, or neither of these. 

(a) a’ © 6.—There is an arc @ D aso that its image @’C 5. By moving 
the points of C; and of C;’ continuously so that @ approaches the whole 
curve C; and a’ a point of C;’, the characteristic remains constant. In the 
limit, J, = 0. 

(b) a’ > 6.—There is an are @ © aso that its image a’ D> 5. Moving the 
points of C; and of C;’ so that @ approaches a point of C; and @’ the whole 
curve C;’, one sees in the limit that J; = —2. 

(c) Neither of (a), (b).—Move the arc a’ continuously on C;’ so that its 
final position is interior to y; this can be done, if necessary, by moving P’ 
across Q or Q’ across P. There isan arc aD asothata’ Cy. Letting a@ 
approach the whole curve C; and @’ a point of y, we see that J; = —1. 

Finally, in (3), there remains the case when g(z) has zeros on C;. Sup- 
pose, e.g., that P’ coincides with P. Then J; = —1, —2 or —2 according 
as Q’ is interior to 4, at.Q, or interior to y. 

Therefore, in all cases for the integral J;, where j > 3, J; < 0. 


k 
Substituting the inequalities derived above in the formula N = P+ oI, 
j=1 


one obtains VN < 2 —1 —1 —1 = —1, contrary to N>0._ It follows that 
f(z) = 2, and the uniqueness theorem is proved. 

Similarly, uniqueness can be established by the above method for other 
types of normalizations and other classes of normal domains. 


1 Courant, Manel and Shiffman, ‘‘A General Theorem on Conformal Mapping of Mul- 
tiply Connected Domains,” these PROCEEDINGS, 26, 503-507 (1940). 

2 After completing the present note, the author learned that this method was indicated 
by Carleman for the case of circular boundaries in a brief note in Comptes rendus, 168, 
843-845 (1919). 

3 C, need only be star-shaped relative to the point 0. 

‘It can be shown that two homothetic convex analytic curves intersect in at most 
two points. If C; were not convex, there may be more than two points of intersection, 
and the conclusion J; S 0 would no longer be valid. For example, if C;,’ intersected Cy 
in 2r points, then I; may have values between (—1—r) and (—1++7), inclusive. 














